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Abstract. We derive an algorithm to determine recursively the 
lap number (minimal number of monotone pieces) of the iterates of 
unimodal maps of an interval with free end-points. The algorithm 
is obtained by the sign analysis of the itineraries of the critical 
point and of the boundary points of the interval map. We apply 
this algorithm to the estimation of the growth number and the 
topological entropy of maps with direct and reverse bifurcations. 

1. Introduction 

Topological entropy, introduced in 1965 by Adler, Konheim and 
McAndrew, [1], is an invariant of topological conjugacy for self-maps 
of an interval. Here, using symbolic dynamical techniques, we derive 
a formula to calculate the lap number of the iterates of the class of 
maps of an interval with free end-points, leading to a straightforward 
estimation of the topological entropy of these maps. The lap numbers 
of this class of maps depends on the symbolic itineraries of the critical 
point and of the two extreme points. 

Let / = [a, b] be a closed interval of the real line M. We consider the 
class J-" of C'^{I) maps f: I ^ I, with a critical point at Xc G {a,b), 
and such that, 

(I) : r(x,)<0, 

(II) : /'(x) > for X < Xc, and f'{x) < for x > Xc- 
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In this paper, the single-humped maps in T are called unimodal 
maps. The maps in the class J-" are not necessarily symmetric around 
the critical point, and the iterates of the end-points of the interval / 
need not to converge to the same limit set. 

The chain rule of differentiation applied to the nth iterate of /, writ- 
ten /" (/° is the identity map), shows trivially that along with /, all 
iterates are also twice differentiable on I. In particular, 

(1) r (x) = f'{r-\x))f'ir-'{x)) ■ ■ ■ nx) 

imphes that critical point of for every n > 1. From (1) and 

condition (II) we have: 

Lemma 1.1. // / G J-", then the critical points of , n > 1, are the 

points X G (a, b) such that /*(x) = Xc, for some < i < n—1. Moreover, 
all the critical points of with n > 1, are maxima or minima, but 
not inflection points. 

It follows from Lemma LI, that the critical points of f^ with n > 1, 
are the pre-images of the critical point Xc up to order n — 1. If f{xc) > 
Xc, f{a) < Xc and f{b) < Xc, then the iterated maps with n > 2, 
have more than one critical point. On the other hand, if f{xc) < Xc, or 
/(xc) > Xc, f{a) > Xc and f{b) > Xc, then the only critical point of 
with n > 2, is Xc. 

A well-known example of a family of maps in J-" is provided by the 
one-parameter quadratic (or logistic) family f^{x) = 4/ix(l — x), where 
/ = [0, 1] and the parameter /i G (0, 1]. For this family of maps, the 
critical point Xc = 1/2 is /i-independent. If /i G (0, 1/2], then the only 
critical point of with n>l, isx = Xc = l/2. 

Two self- maps fi and /2 of the intervals Ii and I2, respectively, are 
topologically conjugate if there exists an homeomorphism h : Ii ^ I2 
such that /2 = h o fi o h~^. The topological entropy of a piecewise 
continuous interval map / is an invariant of topological conjugacy, [1], 
and its topological entropy is calculated through the minimal number 
of monotone pieces or laps of the iterates or by the number of fixed 
points of /"■, [1, 14, 15]. To be more precise, denoting by the minimal 
number of monotone pieces of and by A^(/") the number of fixed 
points of then the topological entropy of / is given by, [15], 

(2) h{f) = lim sup - log in = lim sup - log A^(/") 

n— >oo n n, 

The number s = hmsup^.^^^ i}/"' is called the growth number of /. In 
the case of the logistic family of maps f^{x) = 4/ix(l — x), we have 
/i(/i) = log2 and h{f^) = 0, for fi G [0, 1/2]. 
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An important characteristic of a unimodal map / with positive topo- 
logical entropy is the existence of a semiconjugacy to the tent map with 
slopes ±e^^-^\ This is the Milnor and Thurston classification theorem, 
[14, 11]. The variation of the topological entropy as a parameter of a 
family of maps is changed, implies that the number of periodic points of 
the family of maps also changes, showing the existence of bifurcations 
in the dynamics described by interval maps. This is in general associ- 
ated with a change in complexity in the sense of Li and Yorke's chaos 
scenario [12]. The topological entropy for families of interval maps is a 
way to analyze bifurcations and is also a quantitative evaluation of the 
complexity of the dynamics. 

Several authors have proposed approximating algorithms to calculate 
the topological entropy of unimodal maps with fixed boundary points. 
In the case of the logistic map, techniques based on Markov processes 
have been introduced, [2, 10]. In this approach the spectral proper- 
ties of a transfer matrix with increasing dimension are analyzed. This 
transfer matrix is defined by the images by the logistic map [2, 10]. 
Another technique consists in approximating unimodal maps by piece- 
wise monotone maps and then calculating the topological entropy of 
the piecewise approximation to the map, [3]. All these methods con- 
verge to the topological entropy, however they are indirect and do not 
give information about the topological characteristics of the map, as, 
for example, the lap numbers and the number of critical points. 

A direct method of estimating and calculating the topological entropy 
of an interval map is with the lap numbers of the iterates of a map. 
In [5, 6] a recursive formula to calculate was proposed for the class 
of unimodal maps with fixed boundary points. The proof for unimodal 
maps with a local minimum and fixed boundary points was given in 
[7] . Here, we generalize this result for the larger class of maps J-". The 
main result of this paper is Theorem 4.2, where we derive a general 
formula to calculate recursively the lap number in of the iterates of 
/ G J-". This formula depends on the orbits (itineraries) of the critical 
point and of the boundary points of the map. On the other hand, this 
formula can also be efficiently used to estimate the topological entropy 
and the growth number for unimodal maps obtained by time series 
analysis. 

The proof of Theorem 4.2 is based on the kneading calculus of Milnor- 
Thurston [14] and on a symbolic sequence introduced in [5] which is 
derived from the signed itinerary of the critical point of a unimodal 
map. This new symbolic sequence will be called the Min-Max sequence 
(MMS) of the map. The MMS gives a procedure to determine on the 
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interval / the sequence of maxima and minima of any iterated map 
with n > 1. 

To fix the notation, next, we describe the symbohc approach due to 
Metropohs et al, [13], and recall some basics related to the kneading 
calculus. The Min-Max sequence of a unimodal map and its most 
important properties will be analyzed in the next section. 

Consider a map / G J-". Following Metropolis et al, [13], we assign 
the symbols Jq and Ji to the intervals [a,Xc) and {xc,b], respectively, 
and the symbol C to the critical point Xc- We then associate to every 
{f,x) G J-" X /, a sequence 6/(x) = (0/,n(a;))n>o with entries 6/,„(a;) 
belonging to the alphabet A = {Iq, C, Ji}. The sequence Qf{x), called 
the itinerary of x under /, is defined as follows: 

r Jo if r{x) < x„ 

0/,n(a:) = < C if r(x) = x„ 
[ h if rix) > X,. 

The space of all the itineraries will be denoted by A^'\ where Nq = 
{0, 1, ...}. Since, 

(3) 0/,n(/(x)) = <df,n+lix), 

the action of / on / translates into a left shift on The sequence 

7/ = (7/,n)neN, defined as, 

(4) 7f,n = Qf,ni^c) = Qf,n-lifiXc)), U > 1, 

is called the kneading sequence (KS for short) of /, [14]. Hence, 7/ is 
the itinerary of the critical value f{xc)- 

The signed itinerary of x G / under the iteration of /, is the sequence 
9j-(x) = {Q'f^ni^))n>o on the extended alphabet G {— /i, —C, —Iq, 0, /q, C, Ji}, 
where, 

(5) e^^,o(x) = 0/,o(x), e^,i(x) = e;,i(x), 

(6) Ql„{x) = £(e;,i(x))...£(ey,„_i(x))e/,„(x) for n>2, 
and, 

e(/o) = 1,£(C) = 0, and e(Ji) = -1. 

Note that, e{Io), e{C), and coincide with the signs of /' on Iq, C, 
and Ji, respectively. 

We say that a sequence a G is admissible if there exists f E 
such that a = G^(/(xc)). From now on, we shall always consider 
admissible sequences without explicitly stating it. 

In the following, we write Q{x) instead of Qf{x) whenever the map / 
is understood from the context, and the same applies to the kneading 
sequence and signed itineraries. 
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2. Geometry of the signed itineraries: The Min-Max 

sequenge 

We now relate the signed itineraries of the critical point of a unimodal 
map / G J-" with the structure of maxima and minima of the iterates 
of/. 

Suppose that f^, n > 1, has a local maximum (resp. minimum) at 
some point x ^ I = [a,b]. In order to simplify the notation and the 
proofs, in the following, we say that /"(x) is a "positive" maximum 
(resp. minimum), if /"(x) — Xc > 0. If, otherwise, /"(x) — Xc < 0, 
then we say that /"(x) is a "negative" maximum (resp. minimum). In 
the remaining case, /"(x) — Xc = 0, and we say that /"(x) is a "zero" 
maximum (resp. minimum). 

Lemma 2.1 ([5]). Let f eJ', andk> 1. Then: 

(a) : // f^{x) = Xc and f{xc) > Xc, then f^^^{x) is a positive 
maximum. 

(b) : If f^{x) is a negative [resp. positive] minimum, then f'^^^{x) 
is a minimum [resp. maximum] . 

(c) : If f^{x) is a negative [resp. positive] maximum, then /'^^^(x) 
is a maximum [resp. minimum]. 

(d) : If f{xc) < Xc, then f^{x) has only one critical point, a neg- 
ative or a zero maximum. 

Proof. For > 1 we have, 

(7) f'+''{x) = if o f'nx) = nfix))f"{x) 

and, 

(8) r+i"(x) = r(/^(x))(r (x))2 + f{fH^))f"ix). 

(a) If f'{x) = Xc, by (II), f^+^'ix) = f'{xc)f'"{x) = 0, and by (I), 
fk+^"{x) = nxc){f"{x)f < 0. As f{Xc) > Xc, f'+\x) = f{Xc) > Xc. 

(b) If /'^(x) is a minimum, then f^^^'{x) = /'(/'^(x))/^'(x) = 0, and, 

f'^"\x) = nf\x))f'"{x), 

where f^"{x) > 0. Hence, f^^^"{x) has the same sign as /'(/^(x)). If 
/'^(x) < Xc, then, by (II), f^~^^"{x) > and /'^^^(x) is a minimum. 
Likewise, if /''(x) > Xc, then, by (II), /^+^"(x) < and /^+-^(x) is a 
maximum. 

The proof of (c) is similar and (d) follows from Lemma l.L □ 
Let / G J-" and set, 

(9) iS*^ = {x G / : X is a critical point for f''}. 
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In particular, = {xc}- According to Lemma 1.1, for A; > 1, iS*^ 
contains Xc and its preimages by the iteration of / up to order k — 1. 
This same lemma or Eq. (7) imply that if x G (a, b) is a critical point 
of f^, k > 1, then x is also a critical point of /" for n > k. Hence, 
^ j^fc+i. These critical points can be maxima or minima, and the 
corresponding critical values can be greater than, equal to, or smaller 
than the critical point Xc- 

In order to distinguish all these possibilities, we introduce the new 
alphabet 

(10) M = {m-,m^,m+,M-,M'^,M^}, 

where "m" stands for minimum and "M" stands for maximum. The 
superscript signs attached to m and M specify additionally whether the 
extreme in question is positive, negative or zero in the sense explained 
above in the beginning of this section. 

Now, we define the sequence Uf = {ujf^n)n>i G M.^ as follows: 



m~ 


if /" 




is 


a 


"negative" minimum. 


■m~^ 


if /" 


'fc) 


is 


a 


"positive" minimum. 


mP 


if /" 




is 


a 


"zero" minimum. 


MO 


if 


yXc) 


is 


a 


"zero" maximum. 


M- 


if r 


^Xc) 


is 


a 


"negative" maximum 


M+ 


if r 


yXc) 


is 


a 


"positive" maximum. 



The sequence is called the Min-Max sequence of / G J-", or MMS for 
short, [5, 7]. For example, if /"(x) is a positive maximum, then we say 
that the extremum /"(x) is of type M"*". As usual, we say that coi and 
Uj have opposite signs if the product of their signs is strictly negative. 

The geometric meaning of the MMS of / is clear: The nth iterate of 
/, has a critical point of type G at x = Xc- 

Furthermore, if / G J-", according to (I), 

(11) sign f"{x,) = -1. 

If n > 2 and /^(xc) 7^ Xc for 1 < k < n — 1, then, by (1) and (II) we 
have 

r"{x,) = f{r~\x,)) ■ ■ ■ /'(/(x,))r (xj ^ o, 

and hence 

(12) sign r"(xe) = -£(7/,i) ■ ■ ■ ^(7/,n-i), {n>2). 

Comparison of (11)-(12) with (5)-(6) leads to the conclusion that the 
symbols 6j„(xc) can be identified with the symbols cof^n as in Table 1, 
provided 6j_„(xc) 7^ 0, for n > 1. More specifically, by the rules of 
Lemma 2.1, the symbols Jq, C and Ji in Bj„(xc) translate into the 
superscript signs — , and +, in w/,^, respectively, while the signs + 
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and — in Oj„(a;c) translate into the symbols M and m in ujf^n- If, for 
some k > 2, •yj^k = C (thus £{'jf,k) = and 0/^^+1 (^c) = 0), then the 
construction of the MMS for n > k + 1 follows the rules of Lemma 2.1. 
These rules are summarized in Table 2. 









h 




M+ 


-h 




m+ 


lo 




M- 


-lo 




m~ 


C 






-c 







Table 1. Correspondence between the symbols of the 
signed itinerary of Xc and the Min-Max sequence (MMS). 



Moreover, Lemma 2.1(a-c) also implies that consecutive symbols in 
the MMS obey the transition diagram in Table 2. 













M+ 


m+,M- 













Table 2. Possible transitions between consecutive sym- 
bols of the MMS. 

In the following, we use the shorthand notation (Q;o...a„)°° for any 
symbolic sequence that repeats indefinitely the block (oq. ..«„). A sym- 
bolic sequence of the form (oq . . . a„(a„+i • • • an+k)°^) is called eventu- 
ally periodic. 

Lemma 2.2. Let 7 and u be the KS and MMS of f ^ J-", respectively. 
If 'y is periodic of period k, then u is also periodic and has period k or 
2k. If 'y is eventually periodic, then u is also eventually periodic. 

Proof. Suppose that the KS 7 is periodic with period fc, 7 = (71. ..7^)°°. 
Without loss of generahty assume that f{xc) > Xc, otherwise all the KS 
have one of the forms 7 = (/o)°° or 7 = (C)°°. Thus, yk+i = 7i = h, 
and Ql_^i{xc) = ±7fe+i = ±/i- From Table 1, it follows that Uk+i = M"*" 
or ujk+i = m~^. Since Ui = M"*", in the first case, it follows that u has 
period k. In the second case, since yk+j = 7j for j > 1, we deduce that 
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f'[f^^^[xc)) and f'{f^{xc)) have the same sign for j = l,...,k, hence, 

by (I), 

/''"-'"{x,) = nf\x,)) . . . f'{f'-\x,))nfix,)) ■ ■ ■ /'(/(a;,))/" (xe) < 

and f'^^~^^{xc) is a maximum. By the periodicity of the KS, 72fc+i = 
7fc+i = /i, which means that P''^^{xc) > Xc, thus ijJ2k+i = M^. 

The second part of the lemma follows readily. □ 

Example 2.3. We consider unimodal maps f E r\ C^{I) with nega- 
tive Schwarzian derivative, 



Unimodal maps with negative Schwarzian derivative have some special 
properties, like possessing at most one stable periodic orbit, [16]. 
It can be shown that 

{iMhr, hioihr, ihiohr and h{ior 

are admissible KS for maps of this class, [9]. Using the rules (5)-(6) 
to construct the respective signed KS and the rules in Table 1, the 
corresponding MMS are: 

{M+m'm+M+m+M-M+m+)°^, M+m' {m+ M+)°^ , {M+m-m+)°° and 

M+(m-)°°. 

Let / G J-", / : [a, 6] — )■ [a, b], and let be the set of critical points of 
(see (9)). Let Xk [resp. y^] be the leftmost [resp. rightmost] critical 
point of f'', i.e., 

(13) Xk = miniS^^, yk = maxS'^, 

for k > 1. The relation between the critical points of and /^^^ in 
the interval {xk,yk) is described by the following lemma. 

Lemma 2.4. Let f E and Zk^i < Zk,2 be two consecutive critical 
points for f^, with k > 2. Then, 

(a) : If f^{zk^i) and f^{zk^2) have opposite signs, then there exists 
one and only one z^+i £ {zk^i-, Zk.2) such that f^^^ has a positive 
maximum at Zk+i- Furthermore, f^~^^{zk+i) = f{xc)- 

(b) : Otherwise, there is no critical point of f'''^^ on (2^,1, -2^,2) ■ 

Proof. We may assume that f{xc) > Xc (otherwise, by Lemma 2. Id), 
f^ has only one critical point), (a) By assumption, /'^ is monotone 
on [zk^i,Zk^2]- Suppose f''{zk^i) < Xc < f^{zk^2)- By the Mean Value 
Theorem, there exists one and only one z G (-Zfc.i; Zk^2) such that f^{z) = 
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Xc- Set z = Zk+i- Hence, by (I) and as /(xc) > Xc, f^^^{z) = f{xc) > 

Xc) 

f'^'\z) = f\f\z))f^\z) = f{x,)f'{z) = 0, 

and, f+^"iz) = f"ixc)if'"iz)f < 0. Therefore, has a positive 

maximum at Zk+i- The case, f^{zk,i) > Xc> f^{zk,2) is dealt similarly, 
(b) This assertion is straightforward. □ 

3. Counting laps 

Before analyzing the general case / G J-" without any particular 
boundary conditions, we consider provisionally the following additional 
boundary condition: 

(III): /*'(a) < Xc and f''{b) < Xc, for every k > 1 

If / G J-" and f{xc) < Xc, then the boundary conditions (III) are 
trivially fulfilled and, as mentioned in the Introduction, Xc is the only 
critical point of for all /c > 1. This means that x^ = = Xc for all 
A; > 1, see (13). 

Interesting dynamics is possible if /(xc) > Xc- In this case, the 
boundary conditions (III) imply that the orbits of the endpoints of 
the interval J, a and 6, are confined within the interval [a, Xc). This 
happens, in particular, for maps anchored at the boundary points of 
the interval [a, 6], with /(a) = f{h) = a, as is the case of the logistic 
map. 

Lemma 3.1. Let / G J-" with f{xc) > Xc, and suppose further that 
condition (III) is fulfilled. Let S'^ be the set of critical points of f^ , and 
let Xk and y^, k > 1, be the leftmost and rightmost critical points of f'^, 
respectively. Then, 

(a) : /'^ has a positive maximum at x^ with (i) Xk < Xk-i, for 
k>2, and (ii) f^{xk) = /(xj. 

(b) : has a positive maximum at y^ with (i) yu > yt-i, for 
k>2, and (ii) f'iyk) = /(xc). 

Proof. As / is unimodal, Xi = yi = Xc- 

(a) Let k = 2. Then, by (I) and (II), /^"(xj = /'(/(x,))/" (x^) > 0, 
so f'^{xc) is a minimum. Furthermore, the equation /^'(x) = /'(/(x))/'(x) = 
has one and only one solution on (a,Xc) = (a,xi). In fact, by (I) and 
(II), /'(x) > for X G (a,Xc), and /'(/(x)) = implies that /(x) = Xc- 
By (III) and the condition /(xc) > Xc, we have /(a) < Xc < f{xc), 
and the equation /(x) = Xc has a solution on (a,Xc). This solution is 
unique due to the monotonicity of / on the interval (a,Xc). Hence, the 
unique critical point of on (a,Xc) is X2, and X2 = miuiS^. 
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We claim now that P{x2) is a positive maximum. Indeed, as f{xc) > 
/(xa) = Xc ^ = f{xc) > Xc , 

and, by (I), 

f"{x2) = nx,){f{x2)y<0. 
We proceed now by induction, supposing that (a) is true for k = 
2,3, ...,n — 1. As by hypothesis, /""^(x^-i) is a positive maximum, 
then /""^(a;„_i) > Xc and P~^"{xn-i) < 0. Hence, 

r'(x„_i) = nr'\x^^^))r-"'{xn-l) > o. 

So, it follows that f"'{xn-i) is a minimum. Furthermore, the equation 
/"'(x) = = has one and only one solution on the 

interval (a, a;„_i). Indeed, (i) f^~^'{x) > on (a,x„_i) because Xn-i 
is the leftmost critical point of f^~^, and f^~^{xn-i) is a maximum; 
(ii) by (III), the induction hypothesis and the monotonicity of /"~^ on 
(a,x„_i), f"'~^{x) = Xc has a unique solution on (a, We call 

Xn the unique critical point of /" on (a,x„_i). Finally, to prove that 
f^{xn) is a positive maximum, by (I) and as f{xc) > Xc, 

f {Xn) Xq =^ / {Xfi) f{Xf^ > Xc, 

and, 

r"{xn) = f"{xc){r''Mf<Q. 

The proof of (b) is similar. □ 

Given the KS of a map / G J-", it is possible to sketch symbolically 
and qualitatively the graph of f^, for any A; > 1. Suppose first that / 
obeys the boundary condition (III). We proceed as follows: 

(A) : Fix k > 1 and from the KS of / G J-", construct the first k 
terms of the MMS, uj = {ui = ,ijJ2, ..■,uJk)- 

(B) : Draw two perpendicular coordinate axes and divide the ver- 
tical axis into k rows, corresponding, top to bottom, to the 
iterates f^, 1 < i < k, Table 3. The horizontal axis repre- 
sents the interval [a,Xc]. Divide the horizontal axis into k + 1 
columns, corresponding, right to left, to the critical points Xc = 
Xi > X2 > ■■■ > Xk of Lemma 3.1 and to the left endpoint a. On 
the leftmost column, above a, enter on each row i the sign of 
f (a)-Xe, for 1 < i < A; (i.e., "+" if f (a) > Xc, "0" if f (a) = Xc, 
and "— " if /*(a) > Xc). On the rightmost (or "first") column, 
above Xc = xi = yi, enter on each row i the element Ui of the 
MMS, for i = 1, k. Due to the boundary conditions (III) as 
assumed in Lemma 3.1, the leftmost column of Table 3 contains 
only minus signs. However, dropping boundary condition (HI), 
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the left most column entries can have all the three signs. This 
construction is shown in Table 3. 
(C): On the second column of Table 3, above enter Ui^i on 
row i, for 2 < i < A;. Continue filling out the remaining columns 
above Xj. In the step j, enter Wj-j+i on row i, ioi j < i < k. In 
this way, we get a triangular matrix, with the symbol along 
the secondary diagonal, and each column being the down shift 
by one entry of the next column to the right. This construction 
shown in Table 3, is based on the geometrical meaning of the 
MMS and on Lemma 3.1(a). In view of Lemma 3.1(b), we can 
extend the horizontal axis to include the interval {xc,b], but, 
due to the boundary condition (III), we do not get additional 
information. Up to this point, we know the structure of the 
local extrema of /* at the special critical points Xi, ...,Xk (and 













UJl 


1 












UJ2 


2 










L02 


UJs 


3 


















UJi 




^k-2 


^k-1 


UJk 
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a 






Xs 


X2 


Xc = xi = yi 





Table 3. Ordering of some of the critical points of the 
iterates /* on the interval [a,Xc]. This ordering is de- 
rived from Lemma 3.1, and the construction (A)-(C). On 
the leftmost column, we show the signs of /*(a) — Xc, de- 
termined in this case by condition (III). Note that the 
iterates of the leftmost and rightmost critical points of 
/* have the same signed itinerary as the critical point of 
the unimodal map /. However, the localization of the 
maxima and minima of the iterated map is not com- 
plete. To finish it, we must incorporate the results of 
Lemma 2.4, as explained in step (D). 

(D): In order to determine the location of the remaining critical 
points of /^,..., /'^ on the interval (a, Xc), we analyze the symbols 
of the MMS in row 2 of Table 3. Consider the (only) two consec- 
utive elements of the MMS on row 2, namely uJi and uj2. If they 
have opposite sign, then by Lemma 2.4(a), p has a positive 
maximum at e (^2, Xi). Thus, from row 3 downwards, a new 
column is inserted between the X2- and the xi-columns. This 
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new column contains the symbols, ijJi,0J2, .■.,0Jk-2 (see Exam- 
ple 3.2 below). If the sign of the two consecutive elements on this 
row are not opposite, we are in the conditions of Lemma 2.4(b) 
and no new column is created. Next, go to rows 3, 4, etc. and 
compare row-wise all pairs of consecutive elements. If two con- 
secutive elements on the same row have opposite signs, then we 
proceed as before and, under application of Lemma 2.4(a), we 
insert a new column in-between, beginning on the next row with 
the entry ui, followed by U2,u}s, ... on the remaining entries of 
this new column. The resulting diagram displays the qualitative 
structure of local maxima and minima of the iterates /*, with 
1 < i < k. According to Lemma 1.1, all critical points of /* are 
pre- images of Xc up to order i — 1. 

We call the "MM-table or the Min-Max table of /" a table like the 
one in Table 4, extended to the whole interval [a, b]. 

Example 3.2. Let ] ^ T obeying the boundary condition (III), and 
suppose that the KS of / is 7 = {l\l\C^°° . By the rules of Table 2, 
the corresponding MMS is w = (M"''m~m°)°°. In order to know the 
qualitative shape of, say, f^, we follow the steps (A)-(D) above. The 
symbolic table of maxima and minima is represented in Table 4. Row 4 
of Table 4 provides the necessary information to draw qualitatively the 
graph of From this construction, we can calculate directly the lap 
numbers of for 1 < /c < 4. From Table 4, and counting also laps on 
the right hand side of the critical point Xc, we obtain, £1 = 2, £2 = 4, 
£3 = 8 and £4 = 14. 

The construction just done provides the tools in order to prove The- 
orem 4.2 of the next section, where we drop the boundary condition 
(III). The existence of a simple zero of the equation /"(x) — Xc = 
on (a, x„_i) [resp. on (?/„_i,6)] depends on the signs of /"(a) — Xc and 
/"(x„_i) - Xc [resp. /"(?/„_i) - Xc and /"(6) - x^. If the signs of, 
say, /"(a) — Xc and /"(x„_i) — Xc are opposite, then, as in the proof of 
Lemma 3.1, it follows that there exists a unique x„ G (a, x„_i) such that 
/"(a^n) = Xc. Due to monotonicity, x„ is a simple zero of /"(x) — Xc. 
Otherwise, there is no such solution and so x„ = x„_i, since x„_i is 
the leftmost critical point of This translates mutatis mutandis to 
the half interval {xc,b). In the following, we consider the general case, 
without assuming the boundary conditions (III). 

Example 3.3. We consider the map 

faA^) = + /3 e J^, with 

a = 2.8 and (3 = —0.1, which will be analyzed in detail in section 5. 
This map has a critical point at x = 0. We suppose further that fa,i3{x) '. 
[-(1 + (3), (1 + /?)] ^ [-(1 + /3), (1 + /?)], and a = -b = -{l + 13). The 
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M+ 


1 






M+ 




m 


2 




M+ 


m~ 


M+ 




3 


- M+ 




M+ M+ 




M+ 


4 


a X4 


X3 


X2 




Xc = xi = yi 





Table 4. Structure of the maxima and minima on the 
interval [a, Xc] of the first four iterates of a unimodal map 
/ obeying the boundary condition (III), with KS 7 = 
and MMS u = (M+m-m°)°°. After having 
constructed the table for the first row {k = 1), if two 
consecutive symbols of the row number 1 have opposite 
signs, counting with the boundary conditions, then there 
is a new column beginning in the second row with the 
first symbol of the MMS sequence. From this symbolic 
representation we can count the number of laps of each 
iterated map. Due to the symmetry of the boundary 
condition (III), this diagrams extends symmetrically with 
respect to the critical point Xc = Xi = yi to the interval 
[xc, h] , and the distribution of maxima and minima is 
mirrored. 

first elements of the MMS of are {M^m~rn^ M~ . . .). The signs 
of the iterates f^^pi.'^) ^'^^ fa,pi^)^ with n > 1, are (— , +, —,+,—,■■■ ). 
With the construction done above in the steps (A)-(D), the MM-table 
of the first five iterates of fa,^ is represented in Table 5. The difference 
from this case to the one presented in Example 3.2 is on the boundary 
points. Observe that X3 = X2, since f^j^ has a positive maximum at 
X2 and fa^pi.'^) > Xc', X4 = X3 because /^^^ has a negative minimum 
at X3 and fapi^a) < Xc- The equality follows from a similar 

argument. 

4. Main results 

Consider a unimodal map / G J-", let in denote the minimal number 
of laps of /"■, and let e„ be the number of local extrema of with 
n > 1. Since /" is continuous and piecewise monotone, the laps are 
separated by critical points, and the relation, 

(14) 4 = e„ + 1 

holds. 

Furthermore, let s„ stand for the number of interior simple zeros 
of f"{x) — Xc, n > 1, i.e., solutions of /"(x) = Xc, with x G {a,b). 
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+ 


M+ 






m 


2 




m~ 




M+ 


m+ 


3 


+ 




M+ 




M- 


4 




M- 


m- M+ 




M+ 


5 


a 


X5 = . . . = X2 






Xc = Xi= Hi 





Table 5. Structure of the maxima and minima on the 
interval [a, Xc] of the first five iterates of the unimodal 
map fa,(5 of Example 3.3. The first elements of the MMS 
are {M~^m~m*M~M~^ . . .). Due to symmetry of Z^,^ 
with respect to Xc = 0, the distribution of maxima and 
minima of /^^ on [xc,b] = [0,0.9] is the mirrored image 
of the maxima and minima of ^ on the interval [a, Xc] = 
[-0.9,0]. 



f^{x) 7^ Xc for < i < n — 1, and f^'{x) 7^ 0. In particular, for 
unimodal maps f & J^, we have £1 = 2, ei = 1 and si G {0, 1, 2}. If in 
addition the boundary condition (III) is fulfilled, then si = 2. 

Lemma 4.1 ([5]). Let f E J^. Let Cn be the number of critical points 
of and let Sn be the number of interior simple zeros of f"'{x) — Xc- 
Then, for n > 2, 

(15) e„ = e„_i + s„_i. 

Proof. Indeed, e„ equals the number of sign changes of /"'. Then, Eq. 

(15) follows from the relation /"'(x) = f'W'\x))f''-^'{x). □ 

From (14) and (15), we have, 

(16) En ^n~l (^n Cn— 1 

for n >2. Moreover, induction on n with (15) yields, 

n—l n—1 

(17) en = ei + = 1 + 5j. 

i=l i=l 

For each iterate of a unimodal map, maxima and minima alternate 
along the x direction. This is easily seen in Table 4. In the simplest 
case, the row k of the MM-table of / G J-" contains only the symbols 
M"*" and m~ . Therefore, the laps between two consecutive extrema 
on the graph of f^ always cross the line y = Xc- Whether this also 
happens in the intervals {a,Xk) and iyk,b) depends, of course, on the 
signs of f^{a) — Xc and f^{xk) — Xc on the left side of the interval J, 
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and of — Xc and f^{yk) — Xc on the right side of /. All the four 
possibilities are encapsulated in the relation, 



where, 



(19) 



Sk = ek + l- ak- (3k 



^ ^ if sign(/''(a) - Xc) ■ sign{f^{xk) - < 
^'^ ' 1 if sign(/^(a) - x^) ■ sign{f{xk) - > 



Pk 



if sign(/^(6) - ■ sign(/'=(?/fc) - xj < 

1 if sign(/^(6) - Xc) ■ sign{f''{yk) - xj > , 



In the construction of the MM-table of / of Example 3.3, we have set 
Xj = Xj_i {i > 2) every time /*~^(a) — Xc and /*~^(xi_i) — Xc have no 
opposite signs (see Table 5). The first time this occurs, the type of the 
leftmost extreme of /* is not ui but U2; if p{a) — Xc and f^{xi)—Xc have 
not opposite signs again, the type of the leftmost extremum of /*"'"^ is 
U3. In general, if f^^^{a) — Xc and /*~^(xj_i) — Xc have not opposite 
signs j times, for 1 < z < A; — 1, then the leftmost extremum f^{xk) is 
of type Ui+j (see Table 5). Therefore, sign(/'^(xfe) — Xc) = sign(a;i+j). 
A similar discussion holds for the rightmost extremum f^{yk)- 

In general, the row k contains any symbol from the whole alphabet 
Ai = {m' ,m'^ , M~ , , M~^}. In this case, if some symbol Ui on 
the row k belongs to the set, 

Mb = {m+,M-,m^,M^} 

then, none of the two entries adjacent to Ui (eventually including f^{a) — 
Xc or f^{h) — Xc) can have a sign opposite to the sign of a;,. Every time 
this happens, there are two laps, one on each side of the extremum cor- 
reponding to the symbol Ui G M.hi which do not cross the line y = Xc- 
If all such symbols belonging to A^;, correspond to extrema other than 
the leftmost and rightmost extremun (so these are of types M+ or m~), 
one gets from (18), 

(20) Sk = Ck + l- Ok- (3k- 2bk 

where bk is the number of symbols from A^;, on row k of the MM-table 
of / (see rows A; = 3, 4 in Table 4 for an example on the half-interval 
[a, Xc]). If, otherwise, f^{xk) (resp. f^iyk)) is of a type included in A^f,, 
then we have to set = (resp. (3k = 0) in (20) since the term 26^ 
already accounts for the absence of a simple zero of /'^(x) — Xc in the 
interval (a,Xfc) (resp. {yk,b)). 
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From the above discussion we conclude that Eq. (20) holds without 
restrictions, provided that, 
(21) 

a/c = 



if sign(/''(a) - Xc) 

if sign(/''(a) - 

1 if sign(/''(a) - 

if sign(/^(6) — Xc) 

if sign(/'^(6) — Xc) 

1 if sign(/''(6) - Xc) 



■ sign(wi+A(fc-i)) < 

■ sign(wi+A(fc-i)) > and uji+x(k-i) e Mb 

■ sign(wi+A(fc-i)) > and Wi+A{fe-i) ^ Mb 

sign(a;i+p(fc_i)) < 

sign(a;i+p(fc_i)) > and uji+p(^k-i) e Mb 
sign(a;i+p(fc_i)) > and uji+p(k-i) ^ Mb 



where sign((X'j) = —1,0,1 depending on whether Ui G {m^,M^}, 
{m°,M°}, or {m~^ , M~^}, respectively. The functions X{k) and p{k) 
are recursively calculated as follows: A(0) = p(0) = 0, and for k > 1, 
(22) 

if sign(/''(a) - Xc) ■ sign(a;i+A(fc-i)) < 



m 



p{k) 



A(A;-1) + 1 if sign(/'^(a) 



■ sign(wi+A(fc-i)) > 







if sign(/'=(6) - Xc) ■ sign(u;i+p(fc_i)) < 



p{k - 1) + 1 if sign(/'=(6) - ■ sign(a;i+p(fc_i)) > . 

We can now derive the main result of this paper. 

Theorem 4.2. Let tu = {ujk)k>i be the MMS of f E J-', and set, 

Jn = {'^ < j < n : Uj e Mb} 

where Mb = {m~^ , M~ ,m'^, M'^} . Let {ak)k>i and (f3k)k>i be the se- 
quences of integers calculated recursively from (21) and (22), then 



(23) 



'2 ^ (C+i-j 
ieJn 



t 



/3. 



where = 1 and n>l. 



Proof. Let Np be the number of columns that begin at row p on the 
MM-table of /, and, as before, let 6„ be the number of symbols from 
the set {m"'", M~, m°, M°} on row n of the MM-table of /. Then, 



(24) 



By Lemma 2.4(a), Np = Sp-i. Upon substitution of this equality into 
(24) and (20), we obtain, 

Syi €yi -|- 1 2 ^ ^ ^n—j /^n 
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From (17) it follows now that, 

n — 1 

(25) Sn = 2 + ^Si-2^ Sn-j - an- Pn- 

Finally, with the condition ii = 2 and substitution of (25) into (16), 
we obtain (23). □ 

For example, we can check the formulas (20)-(23) with the geometric 
information contained in Table 5 (extended to the half- interval [xc,6]). 
Example 3.3. Here, ak = f3k because the map /2.8 -o.i of Example 3.3 
is symmetric around the critical point Xc = 0, Ji = J2 = 0, J3 = {3} 
and J4 = {3,4}. So, by Theorem 4.2, we obtain. 



k 


ak = Pk 


bk 




Sk 


4 


1 








1 


2 


2 


2 


1 





3 


2 


4 


3 


1 


1 


5 


2 


6 


4 





3 


7 


2 


8 


5 





4 


9 


2 


10 



where bk is the number of symbols of the set A^f, of the row number k. 
In the case of maps with the boundary condition (III), then ak = 
= for every k > 1 (see the proof of Lemma 3.1 and Table 4). So, 
from (23), we conclude: 

Corollary 4.3 ([5, 7]). Let cu = {uJk)k>i be the MMS of f e F fulfilling 
the boundary condition (III), and set, 

Jn = < j <n: Uj e Mb}- 

where Mb = M", m°, M°}. Then, 

(26) 4+1 = 24 ~ 2 ^ ] (4+l-j ^ ^n-j)j 

where 4 = 1 o.iT'd n > 1. 

The theorem 4.2 and the corollary 4.3 provide fast algorithms to 
calculate lap numbers, thus allowing to estimate the topological entropy 
of smooth unimodal maps in a simple and efficient way. Note that 
these formulas involve the MMS of the map and the itineraries of the 
endpoints of the interval /. 

From the technical point of view, the conclusions of theorem 4.2 
remain true if we change the smoth conditions (I) a (II), by continuity 
and monotinicity conditions on the left and right sides of the critical 
point. 
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5. Computing the topological entropy 

Theorem 4.2 and Corollary 4.3 are the basic results in order to esti- 
mate the topological entropy of interval maps. 

The first example we consider here is the one parameter family of 
logistic maps, Xn+i = f^{xn) = 4//x„(l — defined on the interval 
/ = [0,1] and with /i G (0,1]. Clearly, G J-", for all the values of 
/i G (0, 1], and the critical point of the family is Xc = 1/2, independently 
of /i. 



n 


Ai = 0.875 


= 0.9196433776 


/I = 0.957 


^ = 1.0 




(/i/o/i/i)°° 










{M+m-m+M+m+M-M+m+)°° 


M+m'{m+M+)°° 


{M+m-m+)°° 




1 


2 


2 


2 


2 


2 


4 


4 


4 


4 


3 


8 


8 


8 


8 


4 


14 


14 


14 


16 


5 


24 


24 


24 


32 


6 


38 


38 


40 


64 


7 


58 


60 


66 


128 



Table 6. Lap number for several parameter values of 
the logistic map f\i,{x) = 4/ix(l — x). For each parameter 
value, we show the corresponding KS and MMS. 



To calculate the lap numbers for several values of the parameter /x, 
we choose /i = 0.875, fi = 0.9196433776, fi = 0.957 and fi = 1. For 
these parameter values, the logistic maps have the KS and the MMS of 
Example 2.3 and, in Table 6, we show the first lap numbers calculated 
from corollary 4.3. 

To calculate the growth numbers and the topological entropies for 
this family of maps, we can use theorem 4.2 and corollary 4.3 to calcu- 
late numerically the lap numbers and then estimating limsup^^". As 
the rate of convergence of limsup^y"^ to its limit s is not known, in 
order to have an estimate of the numerical errors \s — £l/^\, we use the 
following lemma adapted from [7]. 

Lemma 5.1 ([7]). Let u = {uk)k>i be the MMS of f E T fulfilling the 
boundary condition (III). Suppose in addition that the MMS has least 
period k > 1, u = {uJiU2 ■ ■ ■UJk)°°, and f^{xc) ^ Xc- Then, the growth 
number (s) of f is the largest real root of the polynomial, 

k-1 

P{X)=X''-^ 
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where, 

_ j 1 if Ei = m" or M+ 
~ I -1 if Ei = m+ or M~ . 

For the proof, See the Appendix A. 

To test the convergence of {in)^^^ to the growth numbers of the maps 
we take the MMS in Table 6. For the first and third MMS, we use 
lemma 5.1, and we obtain, s = 1 and s = (1 + a/5)/2, respectively. The 
growth numbers of the second and forth MMS listed in Table 6 are, 
s = a/2 and s = 2, [5]. In the four cases, the growth numbers and the 
topological entropies are. 



(27) 





: s = 


1 


, h = 





iMhr 


: s = 


V2 


, h = 


|log2 




: s = 


2 


, h = 


log 




: s = 


2 


, h = 


log 2. 



In Figure 1, for the same parameter values in Table 6, we compare the 
behavior of {in}^^^""^ as a function of n with the growth numbers (27). 
In Table 7, we show the error between the exact and the numerical 
values of the growth number and of the topological entropy for the 
maps of Table 6. In most cases, the error for the estimation of the 
topological entropy is acceptable for n ~ 512. The worst estimate in 
Table 7, correspond to the MMS of a dynamics on the Feigenbaum 
period doubling bifurcation cascade, prior to the appearance of fixed 
points with odd periods ([5]). 





As = \s-il/''\ 


Ah = \h~ (log4)/?2| 




0.034 


0.033 




0.005 


0.004 


ihiohr 


0.003 


0.002 











Table 7. Error between the exact and the numerical 
values of the growth number and of the topological en- 
tropy for the maps of Table 6, estimated with the choice 



n = 512. 

In Figure 2, we show the topological entropy of the logistic map 
as a function of the parameter /i and calculated from (2). We have 
calculated the KS and the MMS up to iterate number n = 512, and we 
have plotted \og{£n)/n (~ h), for n = 32, n = 128 and n = 512, as a 
function of the parameter fi. The lap number has been calculated 
from Corollary 4.3. For /i < fi^o = 0.892486416 . . . (the Feigenbaum 
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n 

Figure 1. Behavior of the growth number {in)^^^"^^ as a 
function of n, for several parameter values of the logis- 
tic map f^{x) = 4/ia;(l — x). The parameter values are 
the same as in Table 6. The thin lines show the exact 
values of the growth numbers, calculated in (27). For 
n = 512, the errors in the growth number estimates are 
below As = 0.005, except for the first case of table Ta- 
ble 6, where As = 0.034. 




point), the topological entropy of the logistic map is zero ([6, 7]), and, 
for the case n = 512, this condition has been introduced by hand in 
the plot of Figure 2. Near /loo the convergence of \og{in)/n to the 
topological entropy is very slow. Table 7. This can be seen in the 
approximations to the topological entropy with n = 32 and n = 128. 
For these cases, we have not imposed the condition h{f^) = 0, for 

H < /ioo- 

The monotonicity of the topological entropy for the logistic map, 
[8] and [17], shows the existence of bifurcations as the parameter /i is 
varied, and it measures the complexity of the dynamics of the family of 
maps. In the chaotic region, where the topological entropy is positive, 
these bifurcations are associated with bifurcations of fixed points with 
large periods, which are difficult to detect numerically. The variation 
of the lap number and of the topological entropy as the parameter fi 
changes, shows the existence of bifurcations. In [6] and [7], it has been 
shown that for the logistic family of maps the topological entropy in 
the chaotic region is well approximated by the function. 



(28) hifx) 




if yu e [fioo,r]2] 



2 



c(l-/i)-''-2 



if jje [r]2, 1] 
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where a = 1.82957, b = log 2/ log 5 = 0.4498069, c = 7i/2 = 1.534780, 
d = 1/2, /xoo = 0.892486416, 772 = 0.919643377 and S = 4.6692 . . . 
is the Feigenbaum constant. In Figure 2, we have also plotted the 
approximation (28) to the topological entropy. 



0.7 




0.90 0.92 0.94 0.96 0.98 1.00 



Figure 2. Topological entropy of the logistic map as a 
function of the parameter fi, calculated from (2), for n = 
32, n = 128 and n = 512. For /i < /loo = 0.892486416 . . ., 
the topological entropy of the logistic family is zero, [7]. 
For n = 512, we have imposed the condition = 0, 

for /i < /ioo. The dotted line is the approximation (28) 
to the topological entropy. 



We consider now the interval map Xn+i = fa,p{xn), where, 
(29) /a,/3(a;)=e-'^'+/3, 

and a ^ and /3 are real parameters. This two-parameter family 
of interval maps qualitatively describes the dynamics of an electrical 
circuit with a nonlinear diode showing chaotic behavior. The map 
(29) shows direct and reverse bifurcations when the parameters are 
monotonically changed, [4]. 

The map (29) has a critical point at Xc = 0. As fa,i3 G C°°(M), and 
/a,/3(0) — ~2a^ < 0, the map obeys condition (I). Condition (II) is 
trivially fulfilled, and without loss of generality, we take a > 0. To 
choose an interval I such that fa,p maps I to itself, we take the non- 
trivial case (3 G (—1,0]. As /q,/3(0) = 1 + (3, and, 

/',/3(0) = /3 + e-"'(i+^)'<l + /3 (>0), 

we can choose / = [—{l + (3),{l + f3)]. Therefore, all the maps fa, 13 '■ 
[-(1 + /3), (1 + /?)] ^ [-(1 + 13), (1 + /?)], with a > and /3 > -1, are 
unimodal. As the iterates of the boundary points a = + and 
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-1.0 -0.8 -0.6 -0.4 -0.2 0.0 

P 



Figure 3. a) Bifurcation diagram of the map (29) as a 
function of the parameter (3 in the interval [—1,0] and 
for a = 2.8. b) Topological entropy of the map (29), 
with a = 2.8 and (3 G (—1,0]. The topological entropy 
has been calculated by (2) with n = 512, and the lap 
number has been calculated by Theorem 4.2. The pa- 
rameter value where topological entropy changes from a 
increasing to a decreasing function of /3 corresponds to 
the reversion of bifurcations. 

b = {1 + (3) can be on both sides of the critical point Xc = 0, we are 
in the conditions of Theorem 4.2. Therefore, in order to estimate the 
topological entropy of the map, we have to calculate numerically the 
iterates of the critical point and the iterates of the boundary points of 
the interval /. 

To better understand the relation between the topological entropy 
and the dynamics generated by the family of maps (29), we have calcu- 
lated an estimate for the topological entropy through the exact value of 
the lap number obtained from Theorem 4.2, and we have also computed 
the bifurcation diagram of the dynamical system Xn+i = fa,i3{xn)- In 
Figure 3, we show the bifurcation diagram and the topological entropy 
of the maps (29), for the parameter value a = 2.8 and /3 G [—1,0]. As it 
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is seen in Figure 3, the region where the topological entropy is increas- 
ing corresponds to an increase in the complexity of the dynamics of the 
maps. This is due to the emergence of several periodic points. When 
the reversal of bifurcations occur, the topological entropy decreases. 

6. Conclusion 

We have derived an algorithm to efficiently calculate the lap number 
and the topological entropy of the class of single-humped unimodal 
maps with free end-points. We have shown that the lap number is 
determined by the itinerary of the critical point and by the itinearies 
of the boundary points (Theorem 4.2). The algorithm is based on the 
geometric interpretation of the kneading sequence of a map, leading 
us to a new symbolic sequence {oJf,n)n£N — the Min-Max sequence — 
that contains qualitative information about the sequence of maxima 
and minima of the iterates /"(x). 

Appendix A. Proof of Lemma 5.1 

Proof of Lemma 5.1. As f E J^, fulfilling the boundary condition (III), 
we have ak = (3k = 0, and by (25), we obtain. 



(30) 



n+fc— 1 n—1 
Sn+k Si 2 Sfi-f ^_jr' Si -\- 2 ^n—j 

1 = 1 j&Jn + k «=1 jGJn 

n+fc-1 

S'i 2 Sn-\-k—j ■ 

i=n ji^J^ 



As u has least period k, ui = Uk+i = M+ or ui = uJk+i = M~, and 
therefore, as in Table 2, Uk G Aib- Introducing this condition in (30) 
and as f^{xc) 7^ Xc, we get. 



(31) 
where, 



n+k-l 

Sn+k 2 Sn+k—j 

i=n+l jeJk-i 
k-1 

5^ ^j^n+k—j 



1 if = m or 
-1 if Ei = m"*" or M~ . 



The recurrence relation (31), can be written in the matrix form, 
(32) 







\Sn+k/ 





1 ■■■ 0\ 
1 ■■■ 



: : ■■■ 1 

V/c-l £k-2 £k-3 • ■ • ^1, 




24 



RUI DILAO AND JOSE AMIGO 



and P{x) is the characteristic polynomial of the matrix in (32). The 
solution of the recurrence relation (32) has the form, 

m 

(33) s„ = ^Qp,(n)Ar 

1=1 

where, Aj are the eigenvalues of the characteristic polynomial P{x), Ci 
are constants and Pi{n) are polynomials in n of degree rrii — 1, where 
rrii is the multiplicity of Aj. By (14) and (15), we have, 

n-l 

(34) 4 = 2 + ^ . 

i=l 

Substitution of (33) into (34), gives the result of the lemma. 

□ 
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